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Introduction 


This  report  presents  a  method  of  simulating  the  dynamic  behavior 
of  pulsed  IMPATT  oscillators  and  amplifiers.  The  problem  is  greatly 
simplified  by  using  the  quasi-static  approximation,  which  assumes  that 
the  quantities  of  interest  are  changing  slowly  with  respect  to  the  RF 
period.  Although  full  device-circuit  interaction  simulations  can  be 
performed  to  study  pulsed  IMPATT  behavior,  such  simulations  are  costly 
and  it  is  difficult  to  incorporate  dynamic  temperature  effects. 

The  basic  theory  of  the  quasi-static  approach  is  presented  first. 
The  application  of  the  theory  to  the  cylindrical  resonant  cavity  power 
combining  circuit  is  then  outlined.  The  method  of  characterizing  the 
IMPATT  diode  structures  is  explained  and  the  overall  device-circuit 
interaction  simulation  procedure  is  outlined.  Simulation  results  using 
a  high- low  GaAs  IMPATT  operating  at  10  GHz  are  presented.  The  procedure 
outlined  here  can  also  be  employed  to  determine  the  properties  of  other 
doping  profiles  and  double-drift  structures.  Such  studies  are  presently 
being  carried  out. 

The  Quasi -Static  Approximation  Model 

In  steady-state  ac  circuit  theory,  the  differential  equation 
is  converted  to  a  phasor  equation  by  replacing  the  time  derivative  by 
the  factor  ju).  For  this  case,  the  amplitudes  and  phases  of  the  circuit 
voltages  and  currents  are  constant  with  time. 

The  quasi-static  approximation  assumes  that  the  amplitudes 
and  phases  of  the  circuit  response  vary  slowly  in  time  and  thus  the 
change  is  negligible  during  the  period  of  the  RF  signals.1  Therefore, 
a  circuit  voltage  would  have  the  form: 
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Differentiating  Equation  1  with  respect  to  time  yields: 
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From  Equation  2  it  is  seen  that  the  time  derivative  in  phasor 
space  using  the  quasi-static  approximation  is  obtained  from  multipli¬ 
cation  by  the  complex  frequency  s  =  a  +  jw,  where 


and 


1  dV(t) 

vTtT  dt 
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+  <M>(t) 
dt 


(3) 

(4) 


Equations  3  and  4  may  be  used  as  follows.  It  is  assumed  that 
a  circuit  response  voltage  is  in  a  particular  state  at  time  t,  that  is 
with  a  given  V(t)  and  ij>(t).  Given  this  particular  state,  there  is 
a  unique  value  of  complex  frequency  s  that  solves  the  circuit  equations 
in  phasor  space.  Once  this  value  of  s  is  obtained,  Equations  3  and  4 
predict  how  the  circuit  voltage  amplitude  and  phase  will  change  with 
time.  A  time  step  advancement  At  is  taken  using  Equations  3  and  4  to 
calculate  V(t  +  At)  and  4>(t  +  At).  For  this  new  amplitude  and  phase 
the  circuit  equations  are  solved  for  a  new  complex  frequency  s  and  the 
process  is  repeated.  In  this  way  the  time  evolution  of  V(t)  and  4>(t) 
in  Equation  1  is  determined. 
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Applicat ion  to  IMPATT  Diodes  in  a_  Cylindrical  Resonant  Cavity  Circuit 

An  IMPATT  circuit  which  is  currently  of  much  interest  is  the 

2 

cylindrical  resonant  cavity.  This  circuit  is  used  to  combine  the 
output  of  several  IMPATT  diodes.  The  circuit  is  illustrated  in  Figure  1, 
'./here  Figure  la  shows  the  coupling  between  one  coaxial  module  and 
the  cavity  and  Figure  lb  shows  how  several  modules  are  combined. 

In  this  analysis,  the  interactions  between  diode  modules  are 

not  included;  rather,  a  single  module  as  in  Figure  la  is  simulated, 

assuming  a  particular  model  for  the  resonant  cavity.  Figure  2  shows 

the  complete  equivalent  circuit  model  attached  to  the  IMPATT  diode 

terminals  at  plane  D.  Rs,  Lg  and  are  the  package-mounting 

parasitics.  For  all  simulations  in  this  report,  the  values  Rg  =  0.25  £2, 

L  =  0.2U  nH  and  C  =  0.32  pF  were  used  (these  L  ,  C  values  are  for 
s  p  s  p 

an  X-band  package  used  by  Microwave  Associates,  Inc.).  The  first 
transmission-line  section,  characterized  by  length  lx  and  impedance 
Zj ,  models  the  quarter-wave  transformer  shown  in  Figure  la.  This 
transformer  matches  the  relatively  high  impedance  at  the  cavity  to 
the  low  impedance  at  the  diode  terminals.  The  next  transmission¬ 
line  section  in  Figure  2  corresponds  to  the  half -wavelength  coaxial 
line  in  Figure  la  between  the  transformer  and  the  resonant  cavity. 

Since  the  frequency  of  operation  wanders  during  the  simulation,  these 
sections  are  not  always  \/k  and  X/2  long;  for  this  reason  the  trans¬ 
mission  lines  are  included  in  the  circuit  model.  If  the  frequency 
were  constant,  a  simple  impedance  inverter  would  adequately  model  the 
coaxial  sections  in  Figure  la. 

In  Figure  la  it  is  seen  that  the  cavity  is  in  series  with  a 
matched  terminating  resistor;  this  resistor  is  R  in  Figure  2  and 
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FIG. 


1 

t 


CYLINDRICAL  RESONANT  CAVITY  COMBINER.  (a)  A  SINGLE 

COAXIAL  CIRCUIT  MODULE  ILLUSTRATING  MAGNETIC  COUPLING 

TO  TM  CAVITY,  (b)  TOP  VIEW  OF  RESONANT  CAVITY  SHOWING 
0  1  0 

MODULE  POSITIONS  AROUND  PERIPHERY.  (RUSSELL2) 
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is  set  to  50  ft  for  all  the  simulations.  The  cavity  itself  is  modeled 
as  a  parallel  RLC  resonant  circuit.  At  resonance,  the  cavity  appears 


as  a  large  real  impedance,  given  by  R  in  Figure  2;  for  the  simulations 

Lj 

in  this  report,  R  =  500  ft  was  used.  L  and  C  are  chosen  to  set  the 
desired  resonant  frequency  of  the  cavity  as  well  as  the  cavity  Q.  For 
these  simulations,  a  constant  value  L  =  0.396  x  10_1°hwas  used  and 
the  capacitance  C  was  varied.  This  sets  the  cavity  Q  at  200  for  opera¬ 
tion  near  10  GHz.  C  must  be  chosen  so  that  the  cavity  resonant  frequency 
is  near  the  desired  operating  frequency,  however,  the  two  need  not  be 
exactly  equal.  In  these  simulations,  it  was  found  advantageous  to  set 
the  cavity  resonant  frequency  at  f  =  9.98  GHz  for  10-GHz  operation. 

If  this  is  done,  the  RF  voltage  can  build  up  at  a  lower  temperature 
and  therefore  sooner  after  the  bias  pulse  is  turned  on.  Varying  the 
capacitance  C  in  Figure  2  corresponds  to  physically  tuning  the  cavity  by 
changing  the  probe  position,  which  is  actually  done  experimentally  to 
optimize  the  performance. 

To  carry  out  the  simulation,  the  program  must  calculate  the 
circuit  input  admittance  at  plane  D  in  Figure  2  for  any  complex  fre¬ 
quency  s,  denoted  by  Yc(s).  Also,  the  derivative  of  this  admittance > 
or  dYc(s)/ds,  must  be  calculated  for  use  in  a  Newton  search  algorithm. 

Method  of  Cha racterizing  the  IMP ATT  Diode 

In  addition  to  the  circuit  admittance  functions  Y  (s)  and 

c 

dYc(s)/ds,  the  program  requires  a  large-signal  characterization  of  the 
IMPATT  diode.  The  finite-difference  program  developed  at  the  Electron 
Physics  Laboratory  and  described  in  a  previous  report3  was  used  to 
calculate  this  information.  The  structure  simulated  in  this  report 
is  a  high-low,  single-drift  GaAs  IMPATT  structure  whose  doping  profile 
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is  shown  in  Figure  3.  The  information  required  is  the  diode  admittance 
and  the  dc  voltage  as  a  function  of  RF  voltage,  dc  current  density, 
and  device  temperature. 

An  important  assumption  is  made  at  this  point.  Since  the  circuit 
Q  is  very  high,  the  frequency  of  operation  will  be  close  to  the  reso¬ 
nant  frequency  of  the  cavity  throughout  the  simulation.  For  this  reason 
it  is  assumed  sufficient  to  characterize  the  diode  only  at  the  steady- 
state  operating  frequency  of  10  GHz.  This  is  a  reasonable  assumption 
since  the  diode  properties  do  not  vary  much  for  frequencies  close  to  the 
operating  frequency. 

Large-signal  simulations  for  the  structure  of  Figure  3  were  carried 
out  at  the  following  temperatures:  300,  ^00,  500,  535  and  575°K.  For 
each  temperature,  simulations  were  carried  out  for  three  current  densi¬ 
ties:  J  =  1000,  1250  and  1500  A/cm2.  This  represents  15  sets  of 
simulations.  Within  each  set,  i.e.,  for  each  combination  of  temperature 
and  J  ,  large-signal  runs  were  performed  over  a  range  of  RF  voltages 
from  small  signal  to  a  point  beyond  the  RF  voltage  for  which  maximum 
efficiency  was  obtained.  Also,  a  diode  model  for  J^c  =  0  was  included 
consisting  of  the  cold  capacitance  in  series  with  a  small  resistance  to 
model  the  undepleted  region  of  the  diode. 

These  data  were  tabulated  and  read  into  the  device-circuit  inter¬ 
action  program.  This  provides  the  program  with  Yn(V  ,J,  ,T)  and 

D  nr  dc 

Vdc^VRF,Jdc,T^  f°r  c^^scre'te  values  of  V  ,  Jdc  and  T.  Since  the 

program  requires  knowledge  of  Yn  and  V.  for  any  values  of  VD„,  J  and 

v  dc  nr  dc 

T,  an  interpolation  scheme  must  be  used  between  tabulated  values. 

Figure  k  shows  the  large-signal  results  obtained  for  the 
high-low  structure  at  J^c  =  1250  A/ cm2  for  two  temperatures: 


T  =  400°K  in  Figure  4a  and  T  =  500°K  in  Figure  4b.  It  is  seen  that 

the  range  of  V__,  for  which  simulations  are  carried  out  varies  signi- 
nr 

ficantly  with  temperature;  this  fact  complicates  the  interpolation 

process.  Imagine,  for  example,  that  at  some  point  in  the  simulation 

the  program  requires  knowledge  of  g(V__  =  30  V,  J.  =  1250  A/cm2, 

T  =  450°K).  One  might  be  inclined  at  first  to  simply  average  the 

values  of  G  at  V  =  30  V  in  the  curves  of  Figure  4.  However,  it 
Kr 

is  seen  in  Figure  4a  that  in  order  to  obtain  a  value  of  G  at  V  = 

Rr 

30  V  the  last  values  simulated  must  be  extrapolated,  and  since  the 

-G(  V  „)  curve  has  a  large  negative  slope  at  the  upper  end  of  the  V 
nr  RF 

range,  the  value  obtained  for  G(T  =  400°K,  V  =  30  V)  will  be  very 

RF 

large  and  positive.  This  value  will  be  much  larger  than  the  relatively 
small  negative  conductance  G(T  =  500°K,  VD_  =  30  V)  obtained  from 
Figure  4b  and  it  will  dominate  the  average.  Therefore,  for  VRF  greater 
than  22.5  V  (the  largest  value  in  Figure  4a),  the  extrapolated  conduc¬ 
tance  curve  in  Figure  4a  will  dominate  until  the  temperature  becomes 
very  close  to  500°K.  However,  this  is  not  the  expected  behavior;  the 
value  of  V  for  which  G(V  ,T)  3  0  is  expected  to  vary  gradually  from 

nr  nr 

V0_,  =  22.5  V  to  *  35  v  as  T  varies  from  400°K  to  500°K.  This 

interpolation  scheme  does  not  give  the  expected  behavior  because  it 
assumes  linear  variation  in  a  region  which  is  highly  nonlinear. 

Figure  5  illustrates  the  interpolation  scheme  used  to  circum¬ 
vent  this  problem.  The  quantity  G(V,T^)  is  desired,  where  results  are 
known  for  T  =  400°K  and  T  =  500°K  and  where  400  =  =  500.  If 

RATIO  denotes  the  quotient  of  RF  voltage  ranges  that  were  simulated 
at  the  known  temperature? , 
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ILLUSTRATION  OF  THE  INTERPOLATION  METHOD  BETWEEN  TWO  TEMPERATURES. 


RATIO 


(5) 


G(V,T.)  =  G(Vj,i»00)  +  TFACT[G(V2,500)  -  GCV^OO)]  ,  (8) 

In  this  manner,  a  much  more  reasonable  variation  of  G(v)  with 

temperature  is  obtained.  [If  V,  >  V  ( 4  00 )  or  V,  >  V  (500)  in 
r  1  max  2  max 

Figure  5,  the  same  equations  are  used  but  extrapolated  values  of  G  are 

used  in  the  final  interpolation  in  Equation  8.]  Variation  of  B(V  ) 

nr 

and  Vdc(VRF)  with  temperature  are  calculated  in  the  same  way  as  for 
0(V  ).  Also,  interpolation  between  different  values  of  Jdc  where 

T  is  a  constant  are  handled  by  the  same  method;  for  that  case  the 
vertical  axis  in  Figure  5  would  represent  instead  of  T. 

The  diode  temperature  T  is  allowed  to  vary  as  the  simulation 
progresses.  The  equation  governing  this  variation  is  assumed  to  be 


dT  +  T  ~  Ta  =  rthpd 

dt  T  T 


(9) 
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where  T  is  the  ambient  temperature  (°K),  R  is  the  diode  thermal 
resistance  (°K/W),  is  the  average  dissipated  power  in  the  device 
(w),  and  x  is  the  thermal  time  constant  (s).  For  the  X-band  GaAs 
device  studied  in  this  report,  x  was  assumed  to  be  1.5  us  based  on  the 
results  of  Olson.4  The  dissipated  power  is  found  from 


P 


D 


AREA 


l 


dc  dc 


(10) 


where  AREA  is  in  cm2,  J.  is  in  A/ cm2  and  is  in  mhos/cm2. 

dc  D 

If  Pp  is  assumed  constant  during  the  time  interval  t^  to  t? , 
the  temperature  at  the  end  of  the  interval  is  given  by 

-(t  -t  )/t 

T(t2)  =  [T(t,)  -  T,  -  1  '  ♦Ta.RthPd  .  (11) 

The  dc  current  density  in  the  device  is  assumed  to  start  from  zero  at 
the  beginning  of  the  simulation  and  increase  linearly  until  it  reaches 
a  constant  maximum  value.  Similarly,  at  the  end  of  the  simulation  J^c 
decreases  linearly  from  the  maximum  value  to  zero.  In  the  simulations 
presented  in  this  report,  both  the  linear  rise  and  fall  times  were 
chosen  to  be  1  ns. 


The  Overall  Simulation  Procedure 

This  section  describes  the  method  used  to  carry  out  the  entire 
device-circuit  interaction;  references  are  made  to  specific  portions  of 
the  program  that  were  described  in  the  preceding  sections. 

Figure  6a  shows  the  RF  circuit  configuration.  1^  ( )  is  a 
sinusoidal  injection  current  source  used  to  investigate  the  pulsed- 
locked  oscillator-amplifier  case;  it  is  set  to  zero  to  look  at  the 
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free-running  oscillator  case.  Y  (s)  is  the  circuit  admittance  which 

c 

varies  with  the  complex  frequency  s;  this  is  the  input  admittance  of 

the  circuit  of  Figure  2.  Y  (V  ,J  ,T)  is  the  large-signal  diode 

u  Rr  QC 

admittance,  which  was  determined  for  many  combinations  of  V  ,  J 

RF  dc 

and  T  by  running  a  finite-dif ference,  large-signal  program  as  des¬ 
cribed  in  the  previous  section.  The  simulated  values  are  stored  in 
tables  and  intermediate  values  of  Y^  are  found  by  the  interpolation 
procedure  described  with  reference  to  Figure  5. 

Figure  6b  shows  the  bias  circuit  consisting  of  a  current  gen¬ 
erator  with  linear  rise  and  fall  times.  A  more  complicated  bias  cir¬ 
cuit  could  easily  be  included,  however,  this  simple  circuit  was  used 
for  the  preliminary  investigations.  of  Figure  6b  is  the  large- 

signal  diode  dc  voltage;  it  was  obtained,  stored,  and  interpolated  in 
the  same  manner  as  Y^.  In  addition  to  the  parameters  shown  in  Figure 
6,  the  temperature  T  varies  according  to  Equation  9. 

In  Figure  6a,  V  and  <t>  are  the  amplitude  and  phase  of  the  RF 
voltage  across  the  device.  As  previously  explained,  these  quantities 
are  regarded  as  functions  of  time  and  are  assumed  to  vary  slowly  with 
respect  to  the  RF  period.  The  time  variation  of  these  quantities  is  de¬ 
termined  by  solving  Equations  3  and  b.  To  find  the  complex  frequency 
s  =  a  +  ju  to  use  in  Equations  3  and  b,  the  following  equation  for  the 
circuit  of  Figure  6a  must  be  solved: 


Yc<s)  ♦  VVRF-VT) 


(12) 


where  the  only  variable  is  s  and  all  other  parameters  are  considered 

locally  constant.  The  admittance  Y  (s)  of  the  circuit  in  Figure  2 

c 

can  be  expressed  as 
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(13) 


Y 

c 


(s) 


p3(s) 

p"  'Cs )  * 

4 


where  p3  and  are  polynomials  of  the  third  and  fourth  order  in  s, 
respectively.  The  coefficients  in  these  polynomials  also  depend  on  s 
since  they  include  the  factors  tanh  (si^/c)  and  tanh  (sij/c),  where 
and  are  the  transmission-line  lengths  and  c  is  the  velocity  of 
light.  Using  Equation  13,  Equation  12  can  be  rewritten  as 


F(s)  =  P,(s)  +  Gn  +  JBn  - 


_L  e-J* 


P  (s)  =  0  .  ilk) 

4 


To  find  the  s  value  for  which  F{ s )  =  0,  Newton's  iteration  method  is 
used.  The  program  starts  with  an  initial  guess  for  s.  Next,  F(s)  and 
dF(s)/ds  are  calculated;  in  the  dF(s)/ds  calculation  the  s-dependence 
of  the  coefficients  of  the  polynomials  p3  and  P^  must  be  considered  so 
that  the  expression  for  this  derivative  is  quite  involved.  Then,  a 
better  guess  for  s  is  obtained  using  the  formula: 


F(s) 

dF( s )/ds 


This  iteration  process  is  continued  until  the  magnitude  of  F(s)  is 
sufficiently  small  compared  to  the  magnitudes  of  the  individual  terms 
in  Equation  1^. 

For  the  case  with  a  nonzero  injection  current  IT ,  an  initiali- 
zation  procedure  must  be  carried  out  before  starting  the  simulation. 

It  is  assumed  that  between  bias  pulses  J^c  =  0  but  the  injection  signal 
is  present  as  in  Figure  6a.  Therefore,  a  substantial  RF  voltage 
exists  before  the  bias  pulse  is  turned  on,  even  though  the  diode  is 
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passive.  If  the  simulation  were  started  with  zero  RF  voltage,  the 
circuit  of  Figure  6a  would  go  through  a  complicated  transient  contain¬ 
ing  many  natural  frequencies  until  this  RF  voltage  is  reached.  The 


quasi-static  form  for  V  (t)  in  Equation  1  is  not  adequate  to  describe 

Rr 

this  transient .  Therefore,  the  steady-state  value  of  V  and  <t*  in 

RF 

Figure  6a  that  exists  before  the  bias  pulse  is  turned  on  must  be 
determined,  and  the  simulation  must  be  started  with  these  values. 

The  equation  which  must  be  solved  is 


W 


*  w  *  vjdc 


o) 


(16) 


where  <oo  is  the  frequency  of  the  injection  current  source  and  Tq  is 
the  starting  temperature.  The  procedure  is  to  vary  V  until  the 
magnitudes  of  the  left-  and  right-hand  sides  of  Equation  l6  are  equal; 
then  the  phase  <j>  is  set  so  that  the  angle  of  the  right-hand  side  equals 
that  of  the  left-hand  side.  This  procedure  is  carried  out  in  a  program 
called  INITIAL. 

Typically,  all  the  parameters  for  the  circuit  of  Figure  2  are 
known  except  for  Z^ ,  the  characteristic  impedance  of  the  quarter -wave 
transformer  near  the  diode.  Also,  the  large-signal  characterization 
of  the  diode  results  in  current  and  admittance  per  unit  area,  so  that 
the  diode  area  must  be  specified.  The  determination  of  Z and  diode 
AREA  is  simplified  by  noting  that  at  the  design  frequency,  l  and  ^  in 
Figure  2  are  set  to  \/2  and  xA,  respectively;  therefore,  the  two 
transmission  lines  in  series  can  be  replaced  by  a  single  impedance 
inverter  with  constant  K  =  z\  The  problem  then  reduces  to  finding 
the  proper  K  and  AREA  values.  It  is  desired  to  match  the  diode- 
circuit  combination  at  a  particular  frequency  wQ  and  for  a  particular 


large-signal  diode  operating  point,  specified  by  +  jB^  mho/cm2 

and  V  .  The  equation  to  be  solved  is: 
nr 

(Gd  +  JBd)  AREA  +  Yc(j«oo,K)  =  e‘J*  ,  ( 17 ) 

RF 

where  the  circuit  admittance  is  a  function  of  the  as  yet  undetermined 
impedance  inverter  constant  K.  I  and  $  specify  the  desired  injection 

Ll 

vector  at  the  operating  point;  I  is  zero  for  the  free-running  case. 

ii 

Since  Equation  IT  is  a  complex  equation,  it  yields  two  equations  which 
are  sufficient  to  determine  K  and  diode  AREA  for  the  desired  operating 
point.  This  determination  is  carried  out  by  the  program  DETKA. 

Figure  7  presents  a  flow  chart  showing  how  the  quasi-static 
simulation  is  carried  out.  Once  the  circuit  parameters  are  determined 
except  for  in  Figure  2  and  a  diode  operating  point  is  chosen,  the 
program  DETKA  is  run  to  determine  and  the  diode  area  using  Equation 
IT.  Then  if  a  nonzero  injection  signal  is  present,  the  starting  V 

RF 

and  (j>  are  determined  by  solving  Equation  l6  in  the  program  INITIAL.  For 
the  free- running  case,  the  RF  voltage  is  started  at  0.1  V,  which  is  the 
smallest  value  the  program  allows  V  to  assume.  (This  is  essentially 
the  noise  voltage  which  is  assumed  to  exist,  from  which  the  oscillation 
builds  up.)  At  this  point  the  initialization  is  complete  and  the 
QS  program  is  run,  which  carries  out  the  time  evolution  of  the  solu¬ 
tion.  The  linear  turn-on  of  J.  is  started;  the  large-signal  diode 

dc 

admittance  and  dc  voltage  are  interpolated  for  the  new  using  the 
procedure  outlined  in  Figure  5;  the  new  complex  frequency  s  is  deter¬ 
mined  in  subroutine  SOLVE  which  solves  Equation  12  using  a  Newton 
iteration;  after  s  is  known,  V__  and  <p  are  updated  using  Equations 
3  and  4;  after  determining  the  new  VD„,  a  new  diode  dissipated  power 

Hr 
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FIG.  7  FLOW  CHART  OF  QUASI-STATIC  SIMULATION 


19' 


^diss  ca-*-GU-1-ate<i  according  to  Equation  10;  the  temperature  T 
is  updated  from  Equation  11;  and  finally  the  cycle  is  completed  by 


updating  J,  another  time  step.  (Of  course  J.  reaches  a  constant 
dc  dc 

value  after  the  1-nc  rise  time.)  In  practice,  it  is  necessary  to  run 

the  QS  program  at  least  three  times  for  the  entire  simulation  since 

different  time  steps  are  used  during  the  turn-on  period,  the  time  for 

which  is  constant,  and  the  turn-off  period.  Therefore  the  program 

is  written  so  that  it  is  easily  restarted  from  a  previous  solution. 

The  update  of  temperature  T  involves  the  evaluation  of  an 

exponential  as  seen  in  Equation  11;  however,  since  the  time  step  is 

uniform  throughout  the  simulation,  this  exponential  need  only  be 

evaluated  once.  The  solution  of  Equation  3  for  the  V  update  is: 

Rr 


/RF(t2) 


VRF(tl)e 


o(t2-t1  ) 


(18) 


so  that  the  exponential  must  be  evaluated  at  every  time  ^ince 

o  is  changing.  This  is  carried  out  by  the  function  EXPO(x)  in  the 
QS  program,  which  evaluates  the  standard  series  expansion  for  e  and 
retains  the  number  of  terms  in  the  series  sufficient  to  give  1-percent 
accuracy. 


Simulation  Results  Using  a_  Kigh-Low  GaAs  IMP ATT  and  a  10-ys  Bias  Pulse 
In  this  section  simulation  results  are  presented  both  for  the 
free-running  oscillator  case  and  for  the  case  with  an  injection  signal 
present.  For  the  free-running  case,  the  component  values  used  for  the 
circuit  of  Figure  2  as  well  as  the  diode  characteristics  at  the  steady- 
state  operating  point  are  given  in  Table  1.  Z1  and  diode  AREA  were 
determined  in  DETKA  as  was  explained  in  the  previous  section, 
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Table  1 


Diode  and  Circuit  Parameters  for  the  Free-Running  Oscillator  Case 

Circuit  Parameters 


RL  =  500  a 

*2  =  1.5  cm 

l  =  0.21  nH 
s 

C  =  6.1422  pF 

z2  =  50  n 

R  =  0.25  SI 
s 

L  =  0.0396  nH 

tj  =  0.75  cm 

h  -  ° 

R  =  50  n 

Zj  =  33.9  n 

C  =  0.32  pF 

P 

Diode  Parameters  at 

Steady-Qtate  Operating  Point 

Jdc  =  1250  A/cm2  , 

v  - hh -6  7  ■  V  30  7  • 

f  =  10  GHz 

0^  =  -  32.85  mho/cm2 

,  =  210  mho/ cm2 

RTH  =  1T.38°K/W  , 

AREA  =  2.809  x  10-4  cm2  ,  T  = 

500°K 
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Equation  17.  The  parallel  resonant  circuit  modeling  the  cavity  was  tuned 
to  9.98  GHz,  which  allows  the  oscillation  to  build  up  at  a  lower  tempera¬ 
ture  than  if  the  resonance  is  tuned  to  10  GHz,  the  operating  frequency. 

Figure  8  shows  a  plot  of  the  negative  of  the  circuit  admittance  vs. 
frequency,  together  with  the  small-signal  diode  admittance  vs.  tempera¬ 
ture,  at  J^c  =  1250  A/ cm2.  Since,  according  to  Equation  3,  0  must  be 

positive  for  to  build  up  and  since  a  <  0  at  the  beginning  of  the 
hi* 

simulation,  V  cannot  build  up  until  0  =  0,  which  corresponds  to  the 
RF 

intersection  of  the  two  curves  in  Figure  8.  When  this  intersection 
point  is  examined,  it  is  concluded  that  the  r>  'dilation  will  not  build  up 
until  the  diode  has  heated  up  to  almost  h‘jO°K.  Also,  the  frequency  of  the 
RF  voltage  as  the  amplitude  starts  to  build  up  will  be  approximately 
10.27  GHz. 

It  should  be  noted  in  Figure  kb  that  the  magnitude  of  the  diode 

negative  conductance  is  decreasing  at  the  operating  point  of  V„,  =  30  V. 

Rr 

If  an  operating  point  was  chosen  where  the  negative  conductance  magni¬ 
tude  is  increasing  with  V  ,  the  operating  point  would  bo  unstable. 5 
This  type  of  instability,  which  would  occur  even  if  the  diode  temperature 
remained  constant,  will  be  referred  to  as  an  electronic  instability 
in  this  report. 

The  other  type  of  instability  that  can  occur  is  thermal  instability. 
The  criterion  for  thermal  stability  at  an  operating  point  isfi 


THPdiss 


) 


<  1  , 


(19) 


•T 

o 

where  Tq  is  the  operating  point  temperature.  In  this  report  the  thermal 
resistance  is  assumed  to  be  constant,  so  it  may  be  removed  from  the 
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FIG.  8  CIRCUIT  ADMITTANCE  VS.  FREQUENCY  (s  =  ju>)  AND 
SMALL-SIGNAL  DIODE  ADMITTANCE  VS.  TEMPERATURE 
AT  J  =  1250  A/cm*. 


derivative  in  Equation  19.  The  procedure  to  check  thermal  stability  is 


then  to  perturb  the  temperature  to  Tq  ±  AT,  find  the  new  operating  point 
and  thereby  the  new  dissipated  power  P^  at  these  temperatures,  and 
finally  calculate  a  finite  difference  approximation  to  Equation  19. 

This  procedure  was  carried  out  for  the  operating  point  at  Tq  = 
500°K  described  in  Table  1.  If  the  temperature  T  is  slightly  below 
500°K,  the  program  interpolates  between  the  large-signal  diode  character¬ 
izations  stored  at  400°K  and  500°K;  however,  if  T  is  slightly  above 
500°K,  the  program  interpolates  between  values  stored  for  500°K  and  535°K. 
Therefore,  different  results  are  expected  depending  upon  whether  the 
temperature  is  perturbed  above  or  below  500°K.  Perturbing  the  tempera¬ 
ture  from  500°K  to  490°K  yields  the  following  result: 


Pdiss(T  ■  500)  -  Pdlss(T  -  h9!>) 


TH 


10 


=  0.39 


(20) 


so  that  thermal  stability  is  assured  if  the  temperature  remains  below 

500°K.  Calculating  P, .  at  502°K  gave  the  result: 

diss 


PJ.  (T  =  502)  -  PJ.  (T  =  500) 
„  diss  diss 

RTH - 2 


2.01+ 


(21) 


so  that  thermal  runaway  will  occur  if  the  temperature  increases  above 
500°K. 

Figure  9  shows  the  results  for  the  free-running  oscillator 
(I  =0).  The  time  between  pulses  is  many  thermal  time  constants, 
so  that  the  starting  temperature  was  300°K.  The  dc  current  was  assumed 
to  reach  its  maximum  value  at  the  first  time  step;  for  this  case  simula¬ 
ting  the  1-ns  ramp  makes  little  difference  since  o  <  0  at  the  end  of  the 
pulse  rise  time.  Figure  9b  shows  that  the  temperature  reaches 
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FIG.  9  (a)  RF  POWER,  (b)  TEMPERATURE,  (c)  FREQUENCY 

AND  (d)  DC  VOLTAGE  VS.  TIME  FOR  THE  FREE- 
RUNNING  OSCILLATOR  CASE. 


approximately  ^50°K  after  2  ys;  based  on  Figure  8,  the  RF  voltage  is 
expected  to  build  up  at  this  point,  which  is  confirmed  by  Figure  9a. 

(The  time  step  used  was  1  ns,  so  that  although  the  RF  power  appears  to 
turn  on  in  one  step,  there  are  actually  several  time  steps  during  the 
turn  on. )  After  VR?  builds  up,  the  dissipated  power  pd^ss  the  diode 
is  reduced,  since  RF  power  is  being  generated;  this  causes  the  tempera¬ 
ture  to  increase  at  a  slower  rate.  Figure  9c  shows  that,  again  as  ex¬ 
pected  from  Figure  8,  the  frequency  increases  as  V  builds  up  from 

RF 

small  signal.  After  V  turns  on,  the  frequency  is  nearly  constant  at 

nr 

10  GHz.  Figure  9d  shows  that  a  jump  in  the  dc  voltage  occurs  as  V 

nr 

builds  up.  (it  was  often  observed  in  simulating  nonpunch-through 
structures  at  constant  temperature  that  V^c  increases  with  V^,  before 
decreasing  near  the  point  of  maximum  efficiency. )  The  dc  voltage  con¬ 
tinues  to  increase  slightly  as  the  temperature  increases,  since  the 
ionization  rates  are  decreasing  and  therefore  a  higher  field  level  is 
required  to  generate  the  same  dc  current. 

Figure  10  shows  the  turn-off  characteristics  for  this  case.  Note 
that  the  time  scale  is  much  smaller  than  for  Figure  9.  Jjc  starts 

decreasing  at  t  =  0  and  reaches  zero  after  1  ns.  It  is  seen  that  V__ 

Rr 

and  P  decrease  immediately,  since  the  decrease  is  electronically 
dominated  whereas  the  buildup  in  Figure  9  was  thermally  dominated. 

V^c  in  Figure  lOd  decreases  as  J^c  decreases  and  the  frequency  in  Figure 
10c  assumes  a  lower  value  as  the  diode  admittance  rapidly  changes. 

The  temperature  at  the  end  of  the  simulation  in  Figure  9 
was  498. 2°K,  so  that  thermal  instability  was  not  encountered.  To  test 
Equation  19,  the  simulation  was  restarted  at  a  temperature  of  502°K 
and  thermal  runaway  did  occur  as  predicted  by  Equations  19  and  21. 
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RUNNING  CASE. 


Next,  the  case  with  an  injection  signal  present  was  investigated. 


Table  2  gives  the  parameters  for  this  case.  The  injection  current  I 

L 

was  chosen  so  that  the  power  from  the  injection  source  is  approximately 
10  dB  down  from  the  diode  RF  power  at  the  steady-state  operating  point. 
This  input  power  can  be  found  approximately  from: 


(22) 


The  program  DETKA  was  run  to  solve  Equation  IT  for  K  and  AREA. 


Figure  11  shows  a  plot  of  the  negative  of  the  circuit  admittance  and  the 


large-signal  diode  admittance  vs.  near  the  steady-state  operating 

Kr 

point.  The  two  curves  do  not  intersect  at  the  operating  point  as  for  the 


free-running  case ,  but  rather  the  admittances  are  separated  by  the 
"injection  vector"  which  is  the  right-hand  side  of  Equation  17,  shown  as 


an  arrow  in  Figure  11.  This  vector  has  zero  angle  since  4>  was  chosen 

to  be  zero,  and  its  magnitude  is  IT/(V  »AREA)  =  0.16/(29  x  2.512  x  10_l< ) 

Rr 

=  22  mho/cm2.  At  the  start  of  the  simulation  (Jdc  =  0),  the  program 
INITIAL  must  be  run  to  determine  the  initial  steady-state  amplitude  and 
phase  of  the  RF  voltage.  For  this  case  it  was  found  that  VRF  =  3.5^  V 
and  $  =  -1.2  rad. 


Figure  12  shows  simulation  results  during  the  1-ns  rise  time  as 
Jdc  increases  linearly  from  zero  to  1250  A/cm2 .  In  Figure  12a,  VRF 
starts  from  the  value  determined  by  the  INITIAL  program.  Because  a 


substantial  V  „  exists  across  the  diode  at  t  =  0,  the  RF  voltage  can 
nr 

start  building  up  immediately  without  the  delayed  start  effect  seen 
for  the  free-running  case.  Figure  12c  shows  that  the  frequency  is 
greater  than  the  steady-state  operating  value  during  the  turn-on,  which 


is  expected  since  Bp  for  lower  RF  voltages  is  smaller  and  therefore 


Table  2 


Diode  and  Circuit  Parameters  for  the  Injection-Locked  Case 


Circuit  Parameters 


RL  =  500  ft 

fcj  =  1.5  cm 

L  =0.24  nH 
s 

C  =  6.422  pF 

z2  =  50  0 

r  =  0.25  a 

s 

L  =  0.0396  nH 

=  0.75  cm 

I_  =  0.16  A 

1j 

R  =  50  fi 

Z1  =  46.11  0 

O 

II 

C  =  0.32  pF 

r> 

Diode  Parameters  at  Steady-State  Operating  Point 

Jdc  =  1250  A/cm2  ,  Vdc  =  W.3  V  ,  VRfl  =  29  V  ,  f  =  10  GHz 

Gd  =  -  34.93  mho/cm2  ,  =  209.9  mho/ cm2 

Rth  =  19. 57°K/W  ,  AJREA  =  2.512  x  10"4  cm2  ,  T  =  500°K 
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/ 


the  Y^  curve  of  Figure  11  moves  down  into  the  portion  of  -Yc(f) 
corresponding  to  higher  frequencies. 

Figure  13  shows  results  starting  just  as  J,  reaches  its  constant 

-  ~  dc 

maximum  value  and  continuing  for  1  ys.  For  this  portion  of  the  simulation 

it  was  necessary  to  use  At  =  0.1  ns  rather  than  1  ns  as  for  the  free- 

running  case;  therefore,  this  10-ys  simulation  was  performed  in  steps. 

The  gradual  increase  of  V„  in  Figure  13a  is  due  to  the  gradual  temperature 

nr 

rise  shown  in  Figure  13c.  Note  that  because  the  delayed  start  effect 
does  not  occur,  the  temperature  rises  more  slowly  than  for  the  free- 
running  case  in  Figure  9,  since  the  diode  dissipated  power  is  lower. 

The  frequency  rapidly  reaches  the  steady-state  value,  as  shown  in  Figure 
13b. 

Figures  l4a  and  c  show  that  the  RF  voltage  and  temperature  continue 
to  gradually  rise  during  the  1-  to  2-ys  period.  The  frequency  is  very 
close  to  10  GHz  as  is  seen  in  Figure  l4b. 

Figure  15  shows  results  for  the  remainder  of  the  pulse  width,  from 
2  to  10  ys.  For  this  portion  a  time  step  of  At  =  1  ns  could  be  used, 
since  the  parameters  were  not  changing  as  rapidly  with  time  as  for  the 
earlier  portions.  Figure  15a  shows  that  the  RF  voltage  reaches  the 
steady-state  value  late  in  the  pulse  width.  The  frequency  remains  at 
10  GHz,  as  is  seen  in  Figure  15b.  From  Figure  15c  it  is  seen  that  the 
diode  does  not  reach  the  steady-state  temperature  until  near  the  end  of 
the  pulse.  (Since  the  temperature  never  crosses  500°K,  the  thermal 
instability  problem  discussed  previously  is  not  encountered. )  Figure 
15d  shows  that  the  dc  voltage  gradually  increases  to  the  steady-state 
value  of  44.3  V. 
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FIG.  13  (a)  RF  VOLTAGE,  (b)  FREQUENCY  AND 

(c)  TEMPERATURE  FOR  1  ns  AFTER  BIAS 
PULSE  TURN-ON  FOR  THE  INJECTION-LOCKED 
CASE. 


-33- 


FIG.  14  (a)  RF  VOLTAGE,  (b) 

DURING  THE  1-  TO  2- 


LOCKED  CASE 


TOC  M-SD 

<•) 


TOC  fft-fEO 
<b) 


FIG.  15  (a)  RF  VOLTAGE,  (b)  FREQUENCY,  (c)  TEMPERATURE 

AND  (d)  DC  VOLTAGE  FROM  2  TO  10  ps  FOR  THE 
INJECTION-LOCKED  CASE. 


Figure  l6  shows  the  turn-off  characteristics  as  the  bias  current 

falls  to  zero  during  1  ns.  As  for  the  free-running  case,  Vn„  and  P 

nr  nr 

in  Figures  l6a  and  b  decrease  rapidly  due  to  electronic  effects  (the 

temperature  is  essentially  constant  throughout  the  hO-ns  period). 

Again,  the  frequency  in  Figure  l6c  jumps  to  a  lower  value  as  the  dc 

current  decreases.  There  is  a  jump  in  dc  voltage  in  Figure  l6d,  also 

due  to  the  rapid  decrease  of  J ,  . 

r  dc 

Conclus ions  and  Suggestions  for  Further  Investigations 

The  quasi-static  program  is  seen  to  be  a  useful  tool  in  studying 
the  dynamics  of  pulsed  operation.  Once  the  diode  characterization 

is  completed,  the  simulations  are  relatively  inexpensive  (a  10-us  simulation 
costs  under  $10). 

There  are  many  additional  facets  of  pulsed  operation  that  can  be 
investigated  using  this  program.  One  is  the  question  of  how  the  IMPATT 
diode  operates  under  pulsed  conditions  when  the  dc  current  and  pulse 
width  are  varied  for  a  given  diode-circuit  combination.  For  the  simula¬ 
tions  in  this  report,  the  starting  temperature  was  300°K.  However,  as  the 
time  between  pulses  is  reduced,  the  starting  temperature  begins  to  in¬ 
crease.  Use  of  the  quasi-static  program  in  an  iterative  scheme  would 
allow  the  determination  of  a  consistent  starting  temperature. 

Another  important  study  would  be  the  effect  of  different  IMPATT 
doping  profiles  on  pulsed  operation.  For  example,  low-high-low 
single-drift  and  various  double-drift  structures  could  be  used.  There 
is  probably  a  significant  difference  also  between  operation  using  punch- 
through  and  nonpunch-through  devices  in  each  category.  All  of  these 
effects  are  presently  being  studied. 
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FIG.  16  (a)  RF  VOLTAGE,  (b)  RF  POWER,  (c)  FREQUENCY 

AND  (d)  DC  VOLTAGE  AS  J.  DECREASES  TO  ZERO  IN 

QC 

1  ns  FOR  THE  INJECTION-LOCKED  CASE. 
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